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The effect of a Zeeman magnetic field coupled to the spin 
of the electrons on the conducting properties of the disordered 
Hubbard model is studied. Using the Determinant Quantum 
Monte Carlo method, the temperature- and magnetic-field- 
dependent conductivity is calculated, as well as the degree of 
spin polarization. We find that the Zeeman magnetic field 
suppresses the metallic behavior present for certain values of 
interaction- and disorder-strength, and is able to induce a 
metal-insulator transition at a critical field strength. It is 
argued that the qualitative features of magnetoconductance 
in this microscopic model containing both repulsive interac- 
tions and disorder are in agreement with experimental find- 
ings in two-dimensional electron- and hole-gases in semicon- 
ductor structures. 



A hundred years after the Nobel prize was awarded in 
1902 for the discovery of the Zeeman effect and the subse- 
quent explanation by Lorentz, applying a magnetic field 
continues to be a powerful means to elucidate puzzling 
phenomena in nature. One of the most recent examples 
is the interplay of interactions and disorder in electronic 
systems. This field has witnessed a revival of scientific ac- 
tivity after pioneering experiments in low-density silicon 
metal-oxide-semiconductor field-effect transistors (MOS- 
FETs) found clear indications of a metal-insulator tran- 
sition (MIT) in effectively two-dimensional (2D) systems 
[1]. Until then, electrons in a 2D disordered environment 
were thought to always form an insulating phase; this 
mind-set was based on the scaling theory of localization 
for non-interacting electrons, supplemented by perturba- 
tive treatments of weak interactions, as well as studies of 
the limiting case of very strong interactions. The surpris- 
ing phenomena were soon confirmed in other semicon- 
ductor heterostructures, although the interpretation in 
terms of a quantum phase transition remains controver- 
sial, and a wide variety of experimental and theoretical 
approaches were unleashed at the problem [2] . 

Among these approaches is the application of magnetic 
fields. Contrary to the well-known effect of a magnetic 
field in weak-localization theory to disturb interference 
phenomena and hence undo localization and insulating 
behavior, the negative magnetoresistance effect [3], in 
the Si MOSFETs and similar heterostructures, the mag- 
netic field is found to suppress the metallic behavior and 
therefore promote insulating behavior [4-6]. The effect 



is present for all orientations of the magnetic field rel- 
ative to the 2D plane of the electrons. In particular, a 
Zeeman magnetic field, applied parallel to the 2D plane 
of electrons and therefore coupling only to the spin, and 
not the orbital motion of the electrons, has been used ex- 
tensively. This puts into focus the important role played 
by the spin degree of freedom of the electron, and its 
polarization [7-10]. 

In this Letter, we present a numerical study of a mi- 
croscopic model for interacting electrons in a disordered 
environment including the effect of a Zeeman magnetic 
field. The present study extends our earlier work with- 
out a magnetic field [11], in which we found clear indi- 
cations that interactions enhance the conductivity and 
lead to metallic behavior in a temperature range (about 
1/10 of the Fermi energy) similar to that of experiments. 
Later numerical approaches have sometimes [12-15] led 
to different conclusions from ours, but they either treat 
the problem within a Hartree-Fock method, or else use 
diagonalization methods which deal with considerably 
smaller numbers of electrons than can be studied in our 
approach. Very recently, an improved study using the 
same approach as in Ref. [12] confirmed our main finding 
[16]. While the numerical evidence is mixed concerning 
the occurrence of a MIT due to interactions, there is a 
consensus in favor of a Zeeman magnetic field tuned tran- 
sition [14,15,17,18], as we shall describe in more detail 
below. 

The microscopic model that we study is the disordered 
Hubbard model defined by: 
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where c JCT is the annihilation operator for an electron 
at site j with spin a and n JCT = c^ a Cj a is the occupa- 
tion number operator, tij is the nearest neighbor hop- 
ping integral (i.e. Uj = Q if i and j are not neighbor- 
ing sites), U is the on-site repulsion between electrons 
of opposite spin, /x the chemical potential, and Bu the 
Zeeman magnetic field. We consider a square lattice. 
Disorder is introduced by taking the hopping parame- 
ters Uj from a probability distribution P(£y) — 1/A t for 
Uj G [t— At/2, t+At/2], and zero otherwise. A t measures 
(bond) disorder strength [19]. 
We use the Determinant Quantum Monte Carlo 
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(QMC) method, which has been applied extensively to 
the Hubbard model, both with and without disorder 
[11,19-21]. While disorder and interaction can be varied 
in a controlled way and strong interaction is treatable, 
QMC is limited in the size of the lattice, and the sign 
problem restricts the temperatures which can be stud- 
ied. To alleviate the sign problem, we use off-diagonal 
rather than diagonal disorder, and tune the value of (i 
such that density (n) = 0.5 (where the sign problem is 
less severe). Interestingly, the sign problem is also re- 
duced by the presence of disorder [11]. 

The quantity of immediate interest when studying 
possible metal-insulator transitions is the conductivity 
and in particular its T- and Bn -dependence. By the 
fluctuation-dissipation theorem ad c is related to the zero- 
frequency limit of the current-current correlation func- 
tion. A complication of the QMC simulations is that 
the correlation functions are obtained as a function of 
imaginary time. To avoid a numerical analytic contin- 
uation procedure to obtain frequency-dependent quan- 
tities, which would require Monte Carlo data of higher 
accuracy than can be produced in the presence of even 
a tolerable sign problem and the need for disorder aver- 
aging, we employ an approximation to obtain o~dc from 
the wavevector- and imaginary-time-dependent current- 
current correlation function (see e.g. [11], where also tests 
of the approximation are discussed). Another interesting 
quantity to study in conjunction with the magnetocon- 
ductivity is the degree of spin-polarization P of the elec- 
tronic system: P = (nj, — n-f )/ (n^ + n-f ) , where nj, , nj are 
the average spin-densities of the corresponding number 
operators in (1). 

In order to study the effect of the Zeeman magnetic 
field B» on the metallic behavior, we start from the case 
with density (n) = 0.5 and disorder strength A t = 2.0 
for which the model exhibits clear metallic behavior: o^c 
rising when lowering temperature T [11]. Figure 1 shows 
that turning on B\\ reduces the conductivity and sup- 
presses the metallic behavior; at field strength Bu = 0.4, 
Ode appears T- independent (within the error bars), and 
at larger field strengths shows a tendency to decrease 
upon lowering T. We do not expect t7<j c to go to zero, 
as for a real insulator, unless very low T and very large 
lattices (out of reach of our computational approach) are 
employed. Nevertheless, Fig. 1 shows the qualitative fea- 
tures of a magnetic-field-driven metal-insulator transi- 
tion, similar to what is seen in experiment [4-6]. Previ- 
ous numerical approaches using different techniques have 
also produced this effect [14,15,17]. 

In order to ascertain that we are indeed dealing with 
a critical phenomenon and in order to locate the criti- 
cal field strength, we focus on fields close to B\\ = 0.4. 
It is important to note that the effect of Bn is to po- 
larize the electronic system (with our choice in (1), 
is promoted at the expense of n-j-) and therefore a large 
enough Bu will result in electrons with spin down only 
and, because of the nature of the Hubbard interaction, in 
a non-interacting system [22]. Consequently, in the limit 



of large 2D lattices and low temperature, the hopping 
disorder will force the conductivity to vanish. Subtract- 
ing out the non-zero value of <7d c that we obtain at very 
large Bu is then a systematic way to correct for finite size 
and non-zero T. In Figure 2, we show the resulting <5<7d c 
vs. B\\ for our lowest temperatures. A rather abrupt on- 
set appears of 5ad c below B\\ w 0.5, which agrees with 
the field value where the curves of <7d c vs. T change from 
insulating to metallic (Fig. 1). Our data for a 2D sys- 
tem in Fig. 2 are consistent with a linear vanishing of 
<5<7dc as the (quantum) critical point is approached. At 
present, our results, while presenting compelling evidence 
for the transition itself, are clearly not precise enough 
to obtain critical exponents. Interestingly, a transition 
from insulator to metal upon increasing magnetic field, 
i.e. the known negative magnetoresistance effect, occurs 
in an amorphous three-dimensional Gd-Si alloy (showing 
a MIT at zero field), also with a linear vanishing of the 
conductivity [23]. 

In Figure 3, we show the resistivity p (= l/<7dc) as 
a function of B\\ for low T. The crossing point (B\\ — 
0.35 ± 0.10) demarks a critical field strength B c which 
separates fields for which the resistivity decreases when 
lowering temperature (low-field metallic behavior) from 
fields for which p increases upon lowering T (high-field 
insulating behavior) . It is especially noteworthy that the 
critical field strength (which can be roughly estimated 
to lie between 0.3 and 0.5 from Figs. 2 and 3) is clearly 
lower than the field for which full spin-polarization sets 
in. Indeed, in Figure 4, we show how the spin polar- 
ization P, defined above, behaves as a function of B\\ 
at the lowest temperature used: there is no reflection of 
the critical field strength in the behavior of the polariza- 
tion and full spin-polarization only happens for Bu > 1.2. 
This feature of our data is in agreement with recent ex- 
periments performed on 2D electron- and hole-gases in 
GaAs and AlAs [9,10]. Since complete spin-polarization 
is equivalent to a non-interacting system, the separation 
of the two field strengths and the incomplete polariza- 
tion at the MIT present evidence that the Zeeman field 
tuned MIT must be seen as a property of a fully inter- 
acting many-body system, at least in the 2D disordered 
Hubbard model. 

Another interesting feature of Figure 3 is what ap- 
pears to be the saturation of resistivity at a field not 
much higher than B c . Experiments also show this be- 
havior, but only for AlAs, where the saturation is shown 
to coincide with full spin polarization [9] . We argue that 
the on-site nature of the interactions in the Hubbard 
model make the saturation happen at much reduced field 
strength compared to that of complete polarization: at 
our rather low total density the minority spin species will 
effectively be decoupled from the majority spin species 
and both spin species form non-interacting subsystems 
at a field where the minority spin has not disappeared 
completely. Increasing magnetic field further at constant 
total density will then not change the conducting prop- 
erties anymore. 
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The notion of a predictable and straightforward effect 
of B\\ is also concordant with the phenomenon that p(B\\ ) 
behaves qualitatively the same in the metallic and insu- 
lating phases (see e.g. Ref. [2]), and therefore the same 
physical mechanism seems at play in both cases. Our 
results suggest the reduction of the effective interaction 
by spin polarization as a likely candidate for this mech- 
anism. 

In summary, applying a Zeeman magnetic field in the 
2D disordered Hubbard model reduces the effect of the 
Hubbard interaction and is able to bring about a tran- 
sition from a metallic phase to an insulator at a critical 
field strength. We find this critical field is considerably 
less than the field required for full spin polarization, em- 
phasizing that, for the disordered Hubbard model, the 
metal-insulator transition occurs in a region where a con- 
siderable degree of electronic correlation remains. This 
is in good qualitative agreement with experimental ob- 
servations when a magnetic field is applied parallel to a 
2D electron or hole gas in GaAs- and AlAs-based het- 
erostructures [9,10]. For Si MOSFETs, the general phe- 
nomenon of suppression of the metallic behavior is in 
agreement, but the issue of the critical field being smaller 
than a saturating field is less clear [8]. In earlier work, 
we studied the T-dependence of <Jd c for various At with- 
out a -B-field and showed that the Hubbard interaction 
enhances <7d c and leads at low T to metallic behavior 
that can be turned into insulating behavior by sufficiently 
strong disorder. Our present results concerning the effect 
of a magnetic field are consistent with that conclusion: 
the rather strong interactions that caused the conducting 
phase at disorder strength A t = 2.0 (below the critical 
disorder strength of approximately 2.4 above which the 
system is insulating) without _B-ficld are reduced by a 
B-field which is able to drive the system back to its in- 
sulating phase. The latter is also its natural state in the 
absence of interactions. We believe that this consistency 
indicates that the disordered Hubbard model provides 
a coherent, qualitative picture of the phenomena in 2D 
electronic, disordered systems both in the presence and 
absence of a Zeeman magnetic field. 
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FIG. 1. Conductivity Ode (in units of e 2 /K) as a function of 
temperature T for various strengths of Zeeman magnetic field 
£>H . As By increases, a transition from metallic to insulating 
behavior is seen in Ode- Calculations are performed on 8 x 8 
lattices for U/t = 4 at density (n) = 0.5 with disorder strength 
A t = 2.0 (see text); error bars result from averaging over 
typically 16 quenched disorder realizations. B\\ and At are 
given in units of t. 



FIG. 3. Resistivity p as a function of By for various low 
T. The crossing point provides another estimate for the crit- 
ical field strength. Computational details and units are as in 
Figs. 1 and 2; for clarity the error bars have been omitted, 
but can be estimated from Figs. 1 and 2. 



1.0 



0.3 




0.0 0.4 0.8 1.2 

B l 

FIG. 2. Conductivity with value at very high B-field sub- 
tracted, StJdc = CTdc(B||,r)— <7dc(-B|| = 4, T), as a function 
of By for low temperature T. A sharp onset of conductiv- 
ity is seen at a Zeeman field at which the slope of <7dc(T) 
changes sign in Fig. 1. Computational details and units are 
as in Fig. 1; for clarity, only error bars on T = t/6 data are 
shown; those on T = t/8 data are typically slightly larger (cf. 
Fig. 1). 
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FIG. 4. Degree of spin polarization P = (n| — n|)/(ri| +ni) 
(see text) as a function of By for fixed low T — t/8. The polar- 
ization shows little change through the metal-insulator tran- 
sition and is only 0.31 at the estimated critical field strength. 
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